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Abstract: We study ALP-photon-conversion within strong inhomogeneous electromag- 
netic fields as provided by contemporary high-intensity laser systems. We observe that 
probe photons traversing the focal spot of a superposition of Gaussian beams of a single 
high-intensity laser at fundamental and frequency-doubled mode can experience a frequency 
shift due to their intermittent propagation as axion-like-particles. This process is strongly 
peaked for resonant masses on the order of the involved laser frequencies. Purely laser- 
based experiments in optical setups are sensitive to ALPs in the eV mass range and can 
thus complement ALP searches at dipole magnets. 
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1. Introduction 

(Pseudo-)scalar particles or corresponding bound states are often intimately related with 
the realization of global symmetries in particle physics. They can play the role of conden- 
sation channels of symmetry-breaking condensates as well as occur as (pseudo-)Goldstone 
bosons. A particularly prominent example is given by the axion, the pseudo-scalar pseudo- 
Goldstone boson of a conjectured Peccei-Quinn symmetry Q , which so far provides for the 
only viable solution to the strong CP problem. 

Even though the axion or more generally, axion- like-particles (ALPs), generically de- 
velop couplings to photons by means of a dimension-5 operator, giving rise to an effective 
action of the type 

£p = + \d,<pd^^ct> - \m'<P' + Ig^F^uF^^-" , (1.1) 

where, ^ is the axion- like field with mass m and coupling g, direct axion searches are 
inflicted by a presumably extremely small coupling g. Strongest bounds on this coupling for 
a wide mass range exist typically in the sub-eV regime and are provided by astrophysical 
arguments related to stellar cooling or direct solar observation constraining g to lie 
below g < W-^^GeY-\ 

As the relevant stellar axion generation processes involve momentum transfers in the 
keV range, these bounds are somewhat model dependent and do not immediately apply to 
general (pseudo-)scalar particles with couplings (and corresponding form factors) becoming 
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sizable only at much lower momentum transfers Q. This has inspired the realization of 
many laboratory experiments based on optical probing which are specifically sensitive to 
couplings to ALPs at low momentum transfer, typically at \q\ ~ 0{lfieY). Most commonly, 
laboratory searches for ALPs are based on the mixing between axions and photons induced 
by a macroscopic external magnetic field through the trilinear coupling in Eq. ( |1.1D . For 
example, in polarimetric measurements [Q] , one can determine the ellipticity and rotation of 
a polarized light beam induced within an external field and thereby explore the parameter 
space of the ALP mass and coupling. This mechanism has been exploited e.g. at the 
ALP search at PVLAS and is also planned at BMV Further, photon production in 
strongly inhomogeneous fields 0] as well as high-sensitivity interferometric measurements 
Q have recently been suggested as means of ALP detection. 

Alternatively, one can utilize the weak coupling between ALPs and ordinary matter in 
order to shine photons in external fields through light-blocking walls, in so-called "Light- 
shining-through- walls" (LSW) setups 0. This is possible if probe photons are converted 
into ALPs in front of the light-blocking wall and reconverted into photons behind that 
wall, where in typical laboratory searches the conversion processes are induced by strong 
dipole magnets. Such searches are currently performed e.g. by the ALPS [0], LIPSS [11|, 
GammeV [12|, BMV [13| and OSQAR |14| collaborations, also aiming at cosmologically 
relevant scalar fields or hidden gauge bosons 

In both cases, polarimetry as well as LSW, the decisive parameter for the best ob- 
tainable bounds on the ALP mass and coupling is the product of the field strength of the 
external magnetic field B and its spatial extent L being a measure for the optical path 
length. Typically, the dipole magnets which are employed in these setups provide field 
strengths of -B ~ 0{1 — 10)T extending over a length of L ~ 0(1 — 10)m. By use of cavities 
for the probe beam, the interaction region can be extended by a few orders of magnitude, 
depending on the details of the setup. Also, the idea of resonant enhancement in LSW 
setups 1 17] holds out the prospect of considerably enlarging the effective interaction region. 

On the other hand, the highest field strengths which are obtainable nowadays in a 
laboratory are present within the focal spots of high-intensity laser systems. Current Multi- 
Terawatt lasers achieve peak field strengths of i? ~ 0(10^ — 10^)T, however, naturally at 
cost of the spatial extent of these fields, which ranges from L ~ 0(1 — 10)^m. Nevertheless, 
the parameter BL in the laser focus lies within the same ball park as for the dipole searches, 
calling for proposals of ALP search within high-intensity laser-based setups [^]. 

In addition, the achievable laser intensity has gone up by more than six orders of 
magnitude since the invention of chirped pulse amplification [|^], with the prospect that 
the parameter BL within planned facilities such as ELI [20| will considerably exceed the 
equivalent parameter at dipoles within the near future. This makes high-intensity lasers 
a dedicated tool |21| for fundamentals test of QED nonlinearities |22] most prominently 



comprising e.g. the possible detection of vacuum birefringence p3|, 24, 25, Ea, and 



the Schwinger effect |28|. Strikingly, these searches implicitly constitute also a probe for 



physics beyond the Standard model of particle physics [18, 31] 



At first sight, however, the above mentioned optical techniques, namely polarimetry 
and LSW setups for ALP detection seem to be obstructed by the nature of the elec- 
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tromagnetic field configuration witliin high- intensity lasers. High intensities and thus 
high field strengths can only be attained by pulsed lasers with typical pulse lengths of 
r ~ 0(10 — 100)fs and repetition rates of /rep ^ IHz. Thus, the cavity enhancements 
which are used for polarimetric measurements at dipoles are not easily available for purely 
laser-based setups. 

On the other hand, also the insertion of a light-blocking wall in a purely laser-based 
setup is disfavored: In order to avoid damaging of the wall by the high-intensity lasers, the 
two focal spots for the conversion and reconversion processes would in practice be required 
to be separated by ©(cm). As the spatial extent of the focal spots of lasers is by orders 
of magnitude smaller, the generic angular spread of the beam of ALPs released from the 
first focal spot would significantly reduce the number of ALPs that could possibly hit the 
second spot for reconversion. Even if the angular spread could be minimized, a laser-based 
LSW experiment would demand for a temporally very well-synchronized setup. 

In this work, we suggest another mechanism for high-intensity lasers which does neither 
rely on polarimetry nor on light-blocking walls. Since the electromagnetic field provided 
by the high-intensity lasers varies at a scale which can be of the same order of magnitude 
as the wave length of the probe photon, this photon can experience a frequency shift when 
traversing the focal region of the external field owing to the nonlinear term in Eq. (p..lj). 
Consequently, the detection of such frequency shifted photons could thus point towards 
the existence of ALPs. In the following, we compute this effect quantitatively and discuss 
the required setup and specifications of lasers which are necessary for its detection. 

The paper is organized as follows: In Sects. ^ and |3|, we first give the equations of 
motion for the ALPs and probe photons which we reduce to one spatial dimension for 
simplicity where the formalism developed in ||3^ can be applied and discuss the necessary 
parameterization of the high-intensity laser beams. In Sects. ^ and ^ , we compute the 
photon-axion conversion and back-conversion amplitudes in a specific laser configuration 
and discuss the physical reasons why this setup can lead to a frequency shift for the probe 
photons. Finally, in Sect. ^, we will summarize our findings and estimate exclusion bounds 
achievable for the operational high-intensity laser facility at the Institute of Optics and 
Quantum Electronics in Jena (3^ and the planned Exawatt facility ELI [pO|]. 

2. Equations of motion 

As we are interested in the effects of the nonlinear interaction of laser photons, let us give 
the equations of motion for the photon and the axion field that follow from Eq.(l.l): 



Below, we study these equations in a rather general setup, assuming the interaction of three 
independent electromagnetic fields, which we all presume to be provided by high-intensity 
lasers. As it will turn out later, experimentally there is actually just the need for at most 
two sources. 
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We split up the field strength tensors into contributions of a probe beam and 
two external fields Aj and A'^; in addition, we neglect self- interactions of these fields. By 
coupling to Aj, the probe photons can be converted into ALPs, see Eq. (|2.1j). Successively, 
in Eq. ( |2.2| ), these axions can be reconverted into photons a^^^. through a field A'^. For 
simplicity, we specialize to a one-dimensional setup. In addition, since we do not intend 
to focus on polarimetry later on, we assume the incoming probe photons to be polarized 
along the y axis and to propagate along the positive z axis without loss of generality. 

Under these presumptions, employing Coulomb gauge and the metric g = (+,—,—,—), 
Eqs. (|2.lD and (|2.2|) can now be written as 



^d^_d^,+m^)cl>{z,t) = -g[el{z,t)By{z,t) + bUz,t)E^{z,t)\ (2.3) 

id^-d^)al^,{z,t) = -g [Byiz,t)Miz,t)+Eliz,t)dMz,t)] . (2.4) 

For the external fields, we consider different cases of propagation along the ±z axis or 
orthogonal to the z axis. In the first case of parallel propagation the external lasers interact 
through both their electric and magnetic field components. In the second case of orthogonal 
propagation, the external fields can couple only through either their electric or magnetic 
field component. 

In particular we see from Eq. (|2.3|) that the axion amplitude vanishes trivially in a 



setup where the fields e^^ and Ej propagate both along the +z axis, since then efj^ = —bf^ 
and Ej = Bj. By contrast, for counter-propagation of the two beams, the fields on the 
right-hand side of Eq. (^) add up, since then we have Ej = —Bj. 

If the propagation axis of the external field lies orthogonal to the z axis, either the 
magnetic or the electric field component can mediate the conversion process. However, 
the axion amplitude in Eq. (^j^) is invariant under this choice for a linearly polarized 
external beam. On the other hand, the back-conversion from the axions into photons is 
not independent of whether the axion couples to the electric or magnetic field component, 
due to the asymmetric coupling structure in Eq. ( ^.4] ). This is due to the pseudoscalar 
nature of (f). If we had taken to be a scalar field, then the coupling structure in Eq. 
would be given by interchanging B^ and Ef. 

In the following, we solve Eqs. ( ^.3| ) and (|2.4| ) using the retarded Green's functions 
along the lines of ||3^ for the massive and massless differential operators in the equations 
of motion, 

Gl{z,t) = ijo(mVi2-|z|2)0(t-|z|) (2.5) 
G^{z,t) = hit-lzl) , (2.6) 

respectively. 

The solutions to Eqs. ( |2.3D and ( p.4p naturally depend on the details of the kinematic 
setup. In the following, we choose with hindsight a specific setting for which the conversion 
process of photons into axions and vice versa leads ultimately to a frequency shift of the 
probe beam. If this frequency shift exceeds the natural line width of the probe beam, it 
can constitute a measurable signal, possibly indicating the existence of ALPs. 
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3. Parameterization of the fields 



In order to solve Eqs. (|2.3| ) and (2^) we need a parameterization for the electric and 
magnetic fields of the three laser beams. A good model for the spatial inhomogeneities of 
a focused beam is given by Gaussian beams [34|, which are solutions to the paraxial wave 
equation. 

In consistency with our one-dimensional model, we restrict ourselves to the lowest- 
order contribution in the aspect ratio ^0 = 7^) where wq is the waist size and the 
Rayleigh length of the beam^. The waist size of a Gaussian beam is a measure for the 
transversal extent of the beam at the focus, whereas the Rayleigh length parameterizes 
the broadening of the focus along the propagation direction, cf. Fig. They are related 
through the wavelength of the beam as Zr = To maximize the interaction of the 

lasers, their focal spots should have a sizable overlap. Here, we assume the waist size to be 
minimized at the origin x = y = z = for the probe beam as well as for the two external 
fields. 

Note that for Gaussian beams, the existence of the peak external field strength is, of 
course, also limited by a temporal pulse length r. This scale must be larger than the time 
it takes the probe photons to traverse the external fields^: Text ^ -^r^* and Text ^ Wq^^. In 
the present study, we formally work in the limit of infinite pulse length for both external 
field and probe beam: Text, ^"in 00. Our final result will thus be phrased in terms of a 
transition probability for the photons of the probe field. In practice, as real facilities are 
limited in energy, intensity and power, an optimization of the effect under such constraints 
can typically be expected for all pulse parameters of probe and external field being roughly 
of the same order. 

Under these presumptions, the electric and magnetic field for the incoming probe beam 
propagating along the positive z axis reads 



^ini^' = -^fn(^> = = sin ( ujint - h^z + arctan ( — ) ) , (3.1) 



v/1 + (z/4")2 

where Eij^ is the amplitude of the field, Win the frequency, /cin the wave vector and zj." the 
Rayleigh length of the probe beam, as discussed above. From Eq. ( |3.1j ) it can be seen 
that the Rayleigh length not only characterizes the longitudinal extent of the field, but 
also appears in the so-called Gouy phase shift that a focused light beam undergoes, when 



passing through its focus |Q. In addition, it holds in Eq. (|3.lD that = in vacuum; 
nevertheless, we stick to this notational distinction, since it will simplify the discussion of 
energy and momentum conservation later on. Lastly, the probe beam model in Eq. ( |3.1] ) 
as well as the external beams can in general include also a constant phase factor. However, 
since only the relative phase between the three beams is important, we omit such a phase 
factor in the above definition. 



^At higher orders in the aspect ratio, the beam acquires polarization components which are neglected 
in the one-dimensional calculation and loses transversality. 
^Here and in the following c = h — 1. 
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Figure 1: Spatial overlap of three focused laser beams at the coordinate center x = y = z — 0. In 
this picture, the innermost beam gives the Gaussian probe beam embedded in two external beams 
which propagate orthogonally and transversally to it, respectively. The waist size Wg'^' as well as the 
Rayleigh lengths z!^^ and z™* constitute scales for the extent of the beam foci along the propagation 
direction of the probe beam and thus parameterize the drop-off of the electric and magnetic field 
components of the lasers. The parameter 6*0 = 7^ is the aspect ratio for which typically 6*0 <C 1. 



4. Photon- Axion-Conversion 

We can now calculate the photon-axion conversion by solving Eq. ( |2.3| ). With hindsight, 
we choose the external field to propagate orthogonal (_L) to the z axis and discuss the 
implications of other possible settings later. 

Without loss of generality, we choose the external field to couple through its electric 
field component, and plug in its Gaussian beam form 



E^{z,t) = E^sm{uj^t + iP^)e-^^/<^\ (4.1) 
BUz,t) = 0, (4.2) 



where E± is the amplitude, io_i the frequency, Wq the waist size and 'tp± the phase of 
the external beam. Thus, combining the ALP equation of motion, ( |2.3D with the Green's 



function of Eq. (|2.5| ) and the parameterization of the fields, Eqs. ( |3.lD and (4T), we obtain 
for the axion field: 



{z',t') = ^gE,^E^ j 



00 1 



dt"Jo (m^{t' -t"Y -\z' - z"\'^^ e {{t' - t") - \z' - z"\) 

z" 

X sin(a;int" - haz" + arctan(— )) sm{(jj±_t" + ip±) . (4.3) 
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Here we have used primed arguments for the ALP field (j), in order to indicate that these 
are variables over which we still have to integrate in the back-conversion process later on. 

As typical laboratory scales are many orders of magnitude larger than the spatial 
extents of the beams, it is justified to perform the z" integration from — oo to oo. Following 
the constraints for the respective pulse lengths Text and Tin, which were discussed in the 
previous section, we let also the integration over t" extend from — oo to oo for computational 
simplicity. 

We rewrite the sines as sum of exponentials and substitute t" t' — T in Eq. 
yielding 



\t') = — gSinSx / dz" , 



g-iA:in2"gi{a;in+aj^)t'gj(arctan(2;"/4")+'0x) 



j-CO . 

/ dTJQ(m^JW^^ 

J\z'-z"\ ^ 



g-jfcin2"gi(a;in-w^)t'gi(arctan(2"/4'')-'/'x) 



dTJ, 



\z'-z"\ 



b (my/T^ - \z' - z"|2) e-*('^- 



^//|2 ) g-i(i^m+t^±)T'_ 



+ C.C. 



• (4.4) 



The two integrals over T which appear in the above equation evaluate [36] to 



X e 



i6{\ujin i uj±\ — m)sgn(a;in ± uj±) 

"'ax 

-isgn(a;in±a;x)|z'-z"|fc?x i ^("T- ~ l^in ± ^±)|) ^-\z'~z"\ki 

I ,4- 



(4.5) 



where we have abbreviated fc^ = (win ± uj±y — m? and /c^ = m? — {uji^ it wx)^. 

Considering the combined t' and z' dependence of the axion field (f) in Eqs. ( [4.4| ) 



and (4.5), we can already interpret this intermediate result. The first contribution of the 
integral in Eq. ( |4.5| ) encodes the situation where the axion is on shell and the frequency of 
the outgoing axion is equal to the sum or the difference of the frequencies of the interacting 
laser beams. The outgoing axion then propagates with = Win i u!±_ and wave vector 
k^y.. As already suggested by physical intuition, this can only happen if |a;in±c<j^| is larger 
than the mass of the axion m, as encoded by the theta function in front. In addition, the 
axion wave carries transmitted and reflected parts, depending on the sign of {z' — z"). 

The second contribution of the integral in Eq. (4.5) corresponds to the situation where 
the mass of the axion is larger than the sum or the difference of Win and uj±_ , respectively. 
In these situations, the axion production is off shell (os), and the axion wave decays expo- 
nentially fast with a decay constant of fe^. Naturally, this is not the physical situation in 
which we are interested. 
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Thus, specializing to jwin ±a;_L| > m in the following, we find by combining Eqs. (| 
and (O): 



iz',t') = \gE,^E^ 



_J_^i(ujin+ui^)t' ^-isgn{z'-z")ki^z' ^iip±^ 

g-(2"/«,J-)2gi(-fci„+sgn(z'-2")fci)2"giarctan(2"/4") 



dz" 



-(2"/«)(|-)2 i(-A:i„+sgn(2'-z")sgn(a;in-wx)feax)2"p«arctan{2"/4") 



c.c. 



• (4.6) 



In order to evaluate the remaining integrals over z", it is useful to employ the identity 

gi arctan (z/zr) ]^ 



Vl + [z/Zr^ l-iiz/Zr) 



(4.7) 



In this way, the integration over z" in Eq. ( [4 .61) is Gaussian and can easily be performed. 
The remaining integration over S is then most conveniently written in terms of the Error 



function erf(x) = (2/-^/7r) d5exp(— 5^). Eq. (4^) evaluates to 



{z',t') = -7rzl^^gEi^E^ 



— p sin((win + u}±)t' - sgn(z' - z")k^^z' + ^-i)) 







exp A/ct4" + 



sin((win - u}±)t' - sgn(z' - z")sgn(win - u}±)k~-^z' - ij;±)) 

1 I / \ 



sgn(a;in - uj±) 



X ( 1 - erf [ ^ + ^^^"^Q 1 exp ( Afc^z^^^ + ( ^ 



± + 9 



where we have defined 



A/c+ = -fcin + sgn(z' - z")k+^ 

Ak~ = -kin + sgn(z' - z")sgn(a;in - uj±_) 



(4.8) 



(4.9) 
(4.10) 



In summary, Eq. ( |4.8| ) tells us that the axion wave is composed of two partial waves with 
frequencies cj^ = {ujia^uj±), which both have transmitted and reflected parts corresponding 
to sgn(z' — z") = ±1, respectively. 

We find that each of these partial waves in the case of transmission and reflection is 
tied to a corresponding amplitude, which is a combination of an exponential and an error 
function. The basic effect of this factor is that for given beam parameters the partial waves 
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have maximal amplitude for AA;^ ~ if > Wq and for Ak^ ^ if < and decays 
quickly otherwise. The quantitative impact of this damping term, of course, depends on 
the absolute values of and Wq. As we will see later on, for experimentally feasible z^" 
and Wq, it is reasonable to assert the condition AA;^ ~ 0. In particular, we will also find 
an additional damping term for the back-conversion process below. 

As one can conclude from Eqs. ( [4.9D and ( [4.1[1| ), the origin of this damping is con- 
servation of three- momentum of the photon and the ALP: Only if the momentum in the 
conversion process is conserved to a good approximation, the amplitude of the partial wave 
will persist undamped. Let us emphasize that the momentum of the external beam does 
not enter Afc^ at this point, since the external beam propagates transversal to the z axis 
and thus there is no net transfer of momentum in the z direction. This will become impor- 
tant later on. Note also that this damping factor in practice determines the sensitivity to 
the mass of the axion which is contained in fe^. Le. given two laser frequencies, efficient 
conversion into axions can only happen, if the mass of the axion is next to resonance, such 
that the sum of momenta vanishes approximately^: Ak^ ~ 0. Since Ak^ depends linearly 
on the ALP mass, there is exactly one resonant ALP mass for given frequencies Win and 

UJ±. 

Before we proceed with the calculation of the back-conversion of the ALPs into photons, 
let us determine the resonant ALP masses m in the conversion process. In vacuum {k = w), 
we have the requirements 



Akl = -Win + sgn(z' - z")\/(^in+w±)2 -m2 ~ (4.11) 



Akj_ = -Win + sgn(z' - z")sgn(a;in - UJ±)^/ (win - uj±y - m? ~ . (4.12) 
Eq. ( |4.11 ) is solved in the case of transmission (i.e. sgn(2;' — z") = +1) by choosing m = 



w^ + 2winW_|_. For (4.12), there is in principle the resonant solution m = yw^ — 2winW_L 
which implies ^u>± > Win for positive axion masses implying refiection (i.e. sgn(2:' — z") = 
— 1) but requiring also a negative value for w^ = Win— w_l. Thus, the latter is an unphysical 
solution. 

Let us summarize our findings from Eq. ([4.^): induced by the interaction of the probe 
field Win with the external beam wj_, one obtains transmitted and reflected axion waves 
with frequencies Wax = "^in i uj±. For our purposes, we focus on the transmitted partial 
wave with frequency Wax = ^in + ^± since it acquires an undamped amplitude for ALP 
masses m which are close to a resonant mass 



m_L = a/w^ + 2winW_L . (4.13) 



^The above conversion properties are reminiscent of the processes of sum-frequency and difference- 
frequency generation, known from Nonlinear Optics, see e.g. |^^. For these processes, a suitable medium 
with nonlinear dielectric permittivity is used to produce light beams whose frequency is equal to the sum 
or difference of the frequencies of the input beams. There, an efficient conversion can only happen, if so 
called phase-matching conditions are fulfilled. These phase matching conditions are analogous to the above 
condition of three-momentum conservation. 
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Thus, for the following calculation of the reconversion process, we employ for clarity only 
the transmitted axion wave with frequency uj'^^: 



sin((tJin + uj±)t' - k^^z' + 7/;^)) 



zf Ak+w^ 



1 - erf + 




• (4.14) 



Let us finally remark that the damping factor of the amplitude encountered above is in 
fact not an artifact of the Gaussian beam form. This can easily be checked by omitting all 
factors containing Wq and in Eq. (^^), which amounts to calculating the interaction 
between two plain waves. Then, integrating over z" over a length L of the interaction 
region, one finds that the amplitudes of the axion partial waves are proportional to a 
factor of sin(AA;^L/2)/(Afe^), respectively. Hence, also for a plain wave approximation, 
the amplitudes are peaked around Ak^ ~ 0. 

From this ansatz, the conversion amplitude for the ALP in a temporally and spatially 
constant external field follows in the limit uj±,k± — )• 0. The square of this amplitude 
is given in Eq. ( |6.6| ) in the limit^ m ^ cjjn, and will be used later on for a qualitative 
comparison to dipole LSW experiments. 

5. Axion-Photon-Conversion 

We now turn to the back-conversion of the ALPs into photons, by virtue of Eq. (p.4|). 
With the Green's function in Eq. (^) and using Cout = —dtaont, we have to evaluate 

poo pco 

eontiz,t) = -g I dz' I dt'5 {{t - t') -\z- z'\) 



00 J —00 



[[Bl{z',t')d,<P{z\t')]+El{z',t')d,.cP{z',t% . (5.1) 



As in the previous section we choose the external field for the back-conversion with hind- 
sight: Assuming propagation of the electromagnetic wave along the negative z axis, i.e. 
counter-propagating to the transmitted axion wave, we set 



^11 



z 



E^{z,t) = — =^1^^= sin (jyt + /cyz - arctan -[[ + V'H (5-2) 



^Jl+[z/z!) 



E\\ 
1 + {z/z\ 



Bl{z,t) = -^jj(z,t) = ^ ^" = sin I + fcyz - arctan ( ^ ) + ) , (5-3) 



*Note that in this situation for large axion masses m ~ a;in, the computation of the conversion proba- 
bihties requires great care, as discussed in detail in |^^. However, generically, the most stringent exclusion 
bounds on ALPs for constant external fields are obtained for masses m <C itii-a- In the end, we will compare 
the discovery potential for the purely-laser based setup to these bounds. 
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with beam parameter definitions as above. 

In the following it is convenient to use the axion and laser wave again in complex 



notation. Plugging the approximate transmitted axion wave from Eq. (4.14) and the 
external beam (Eqs. ( |5.2D and ( |5.3| )) into Eq. ( ^.ID , evaluating the derivatives acting on 
(j){z,t) and integration over t' , we find 







4 



' yJl + {z'/zlY 



+ C.C. 



l + (2'Ar) 



(5.4) 



Here we have defined 



A/c+ = -/ci + /cy + sgn(z - z'){oj;^ + uj± + uj\\) (5.5) 
A/cij" = — A;^ — k\\ + sgn(2 — z'){uj\a + w_l — f^ii) . (5.6) 



Eq. (5.4) resembles the situation of the ALP production, as the outgoing electromagnetic 
wave Bout essentially consists of two partial waves with frequencies uj^^^ = uj\,^ + uj±_ it uj\\ . 
Again, each partial wave has a transmitted (sgn(z — z') = 1) and a reflected (sgn(2: — z') = 
— 1) contribution. In order to determine the corresponding amplitudes of the partial waves, 
it is necessary to perform the remaining integration over z' . To this end, we make use of 
the first identity in ([4.7D and perform the spatial integration over a closed contour in the 



complex z' plane. In this manner, we find for the integrals in Eq. (5.4) 



J — c 



az — p. \4Je II 

1 + ( z'/z}^ ' 



r dz' J -e'^'^' = vrzl (1 ± sgn(AA;±))e-"" l^'ff I . (5.7) 

■''-oo l±i[z'/zii\ 



Analogous to the previous section, the partial wave amplitudes in Eq. (5.7) entering Cout are 
strongly peaked at vanishing ^k^- As in the photon- ALP conversion, this peak structure 
can be understood in terms of three-momentum conservation, cf. Eqs. ( |5.5| ) and (^.6|), 
where the last term incorporates the momentum of the outgoing photon. 

In addition, sharp cutoffs arise for AA;^ < and A/c^ > in Eq. ( |5.7D through 
the signum function, respectively. Physically, this is due to the Gouy phase anomaly 
1 35]; this behavior is also well known in the context of nonlinear interactions with focused 
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Gaussian beams in media, see e.g. |]37[| . However, one can check by numerical integration 
that the sharp cutoff of the signum function is in fact washed out for integrations over 
finite interaction regions. Thus, for finite, physical interaction regions, the amplitudes are 
maximized for — 0. 

Before W6 continue with, the evalucition of Cout? 

let us check the compatibility of the 
three-momentum conservations for the two conversion processes. Only if the conservation 
of three-momentum is obeyed at both conversions at the same time, the overall amplitude 
is undamped. In addition to conservation of three-momentum, we want to make the im- 
portant additional requirement that Wout / win, as signature for the conversion processes 
to have taken place at all. 



To this end, we consider Eqs. (5.5) and ( ^.6^ again in vacuum (i.e. ui = k). In order 
to compare with the conservation of momentum in the photon-axion conversion process 
encoded in A/c^ (see Eq. ( [4. Ill )), we multiply AA;^^ by —1. This is justified, since the 
exponential damping depends only on the modulus of /S.k^ . Eqs. ( |5.5| ) and ( [5.6[ ) then 
read: 

— AAjjl" = a/ (win -|- ljJi_Y — Tn^ — sgn(z — z')('^m + ^1- + '^ll) ~ "^11 ^ (^•^) 
— AA:|j" = \/ (win + w_l)^ — m? — sgn(2; — z'){uj\-a + lo^ — a;||) -|- ~ . (5.9) 

By comparing the above conditions for A/c^ with that for AA;^ in Eq. ( [4.11D , we notice 
that we need the transmitted part of the outgoing wave, corresponding to sgn(z — z') = 
also for the back-conversion process. Otherwise the required relative sign between and 
the axion wave vector in order to satisfy Eq. ( [4.11 ) can not be reproduced. 



In order to fulfill the condition for A/i;j| in Eq. ( ^.8]) and the condition for AA;J| at 
the same time, a negative frequency contribution of either or would then be needed 
which is clearly unphysical. For this reason we drop this option in the following. 

By contrast, it can be seen that it is possible to obey momentum conservation for 
the conversion and back-conversion processes simultaneously via Eq. (^^): By setting 
uj± = 2a;||, we see that AA;^ = AA;^ = if the axion mass satisfies Eq. (|4.13D . In addition, 
the outgoing electromagnetic wave in this situation has a frequency Wout = '^in + = 
Win + ^\\ which is different from the frequency of the incoming wave. We conclude that also 
our second requirement for Cout is met. 

At this point, it also becomes clear why it is crucial for the frequency shift of the 
outgoing wave to arise that the momentum of the first external laser beam does not enter 
the requirement for /^k~\_ in Eq. ( |4.9| ): Depending on the relative signs of the momenta 
in the two conversion processes, conservation of momentum could only be achieved in 
both processes for io±_ = u)\\. However, this would immediately imply ojia = ^out being 
problematic for experimental observation and require again an LSW setup which is difficult 
to conceive for high- intensity fields. 

In total, it is the above "mismatch" of momentum conservation (as encoded in AA:_|_/||) 
and energy conservation (as encoded in the difference ujont — '^in); which is essential for a 
net transfer of energy to the outgoing photon, favoring experimental detect ability. 
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To summarize, it is the transmitted (T) contribution of the second partial wave in 



Eq. (5.4) for which both respective momentum conservation conditions can be fulfilled 
simultaneously and a frequency shift with respect to the incoming wave arises. Combining 
Eqs. (^.4]) and Eq. (|5.7D and substituting uj± = 2uj\i , we focus for the remainder on 




eZliz,t) « -^g\'zl-4Eir,E^E\\ (^1 - erf + ^^^^ exp (^Akz^ + 

X (1 - sgn (Ak)) e-^"l^'=l ff}l^±^ + A sin{{co,^ + u;\\){t - z) + - ^|;\\) , (5.10) 

where we have set Ak^ = = 

In practice it is of course not directly experimentally assessable which of the two ex- 
ternal beams mediates conversion and which one back-conversion. So far, we have assumed 
conversion to be induced by the field with _L orientation, and back-conversion due to the 
counter-propagating || field. The result for the outgoing wave, however, will depend on the 
order of interaction. As demonstrated in the Appendix, the conversion process with the 
opposite order _L ^\\ results in an outgoing wave with frequency Wout = ^in — ^\\, if one 
chooses UJ± = 2cj||. The resonant mass satisfying momentum conservation in both conver- 
sion processes is then given by my = 2^a;inW|j, being different from the resonant mass at 
interchanged interaction order. In consequence, it is possible to probe the axion coupling 
space around two resonant masses within one setup. 

Our calculations suggest the following experimental setup. Probe photons, which tra- 
verse one counter-propagating and one perpendicularly propagating laser field with fre- 
quencies and uj±, respectively, can experience a frequency shift due to ALP-photon 
mixing. This happens if the frequency of the external perpendicular laser has twice the 
frequency of the external counter-propagating laser: 2a;|| = uj±. The requirement for this 
process is the existence of ALPs with masses close to one of the two resonant masses, which 
are a function of the involved laser frequencies. 

These resonant masses are of the same order of magnitude as the involved frequency 
scales of the lasers, which in optical setups corresponds to ~ 0(eV). As this mass regime 
is so far largely unexplored in laboratory ALP-searches, the proposed experimental setup 
can be complementary to the search involving dipole magnets, see below. 

It is worth emphasizing that the required frequency ratio for the two external lasers 
is fact an enormous experimental advantage, since it implies that indeed only one high- 



intensity laser is needed as external field, since frequency doubling [38| is a standard tech- 
nique even for high-intensity lasers. In addition, the corresponding beam parameters of the 
frequency doubled beam, such as the focal area can in principle be tuned independently 
by the use of appropriate lens systems. 

6. Exclusion limits 

Let us explore the parameter range in the ALP mass and coupling plane which can be 
probed within the presented setup. 
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The number of photons in the beams is proportional to the square of the field am- 
phtudes, being a function of time. However, as the pulse lengths in consideration imply 
a large number of wave trains, a good approximation of the number of frequency shifted 
outgoing photons A'out as a function of the number of incoming photons can be read 
off from a comparison of Eq. ( |5.10| ) (or ( A.15| )) with the incoming field ein in Eq. (p.l|). 



A'out(^^in ±^||) iVin(Win)iVshota 



2 

± • 



(6.1) 



The parameter iVshot counts the number of laser shots used for a measurement. It is 
determined by the total measurement time of data accumulation times the repetition rate 
of the lasers. The quantity is a measure for the probability of the photon-axion-photon 
conversion. Here, a+ denotes the conversion amplitude for outgoing photons of frequency 
Wout = Win + W|| whereas a_ is the conversion amplitude for photons with frequency Wout = 
cjin — uj». They read: 



a+ 



1 



—g\^CzfiE^E\i I 1 - erf ( ^ + 



16 



X exp A/cz™ + 




Win + 2CJ|| 



+ 1 



Uin + 2a;||)2 




a_ = ^^2^3/2^ 

8 4" + zil 



l-sgn(AA;))e-^'-'^*^l , 

(Win +C<J||)2 - m2 



(6.2) 



1 - sgn{5k))e^^''''' + (1 + sgn{5k))e 



-Skz 



e 4 



(6.3) 



where we have inserted the respective axion wave vectors for clarity. The parameters Ak 
and 5k reduce to 



AA; 



6k 



<^in + J (win + 2a;||)2 — 



-Win + Wll + J (cJin + W 



(6.4) 
(6.5) 



As already denoted above, in order to experimentally assess the induced frequency shift, 
it is sufficient to have just one high-intensity laser as external field in this configuration. 
One part of the external beam has to be frequency doubled, such that uj± = 2a;||, while a 
delay line must ensure the simultaneous overlap of the focal spots of the fundamental and 
the frequency doubled beam mode as well as the focal spot of the probe beam. 

In addition, we demand that the three frequency modes satisfy Wjn + wy oo in 
order to make the frequency shifted photons at ujont distinguishable from the frequency 
components already used in the process. We will detail on this requirement below. 

Further, it appears advisable in practice to realize the setup with a slight deviation 
from the exact orthogonal geometry in order to facilitate the detection of the frequency 
shifted photons off the main optical axis in order to reduce the noise amplitude^. 



'For a quantitative discussion of a non-orthogonal setup, the above calculation, of course, has to be 
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Let us now determine the discovery potential in the ALP mass-couphng plane for an 
operational high-intensity facility. First, we discuss a possible setup at Jena |3^, |3^. In the 
near future, the Multi-Terawatt class laser JETI and the Petawatt class laser POLARIS 
pO| ] can be focused simultaneously and synchronized into a single target chamber. Thus we 
can employ the lower-intensity laser JETI for providing the probe photons, while POLARIS 
is used to create the two external fields for the conversion processes. 

Let us first consider the focal parameters of the two lasers. 

To achieve maximum field strengths for POLARIS and a good bunching of the probe 
photons of JETI, we need very small focal spots on the order of the diffraction limit. 
To be more precise, if we define the effective diameter of the focal spot to contain 86% 
(corresponding to the l/e'^-criterion) of the focused beam energy, see below, one has the 
estimate wq ~ /"^A, with the so-called /-number of the focusing lens, which 
characterizes the ratio of the focal length and the focusing aperture diameter. Ambitious, 
but feasible values for the /-number can be as low as /* = 1. Recall that the corresponding 
Rayleigh length is obtained from the waist size through Zr = -j^. 

The laser system POLARIS is designed to provide a peak power of around P = IPW 
(attained through 150J at pulse lengths Text — ISOfs), optimized for a central wavelength 
of Aext = 1035nm, corresponding to = 1.20eV and thus uj_i = 2.40eV. In consequence, 
we find Wq ~ 5^y^ ~ 13.1eV~^, for the frequency doubled beam where we have chosen an 
/-number /# = 5 in order to fit the probe photons into the external field, see below. 

Further, we estimate the Rayleigh length for the counter-propagating fundamental 
beam as zf- ~ vr(/'^)'^Aext — 16.4eV^^ for an /-number of /* = 1. 

The obtained intensities within the focal spot yield /y = 0.86^^ = 2.05 x lO^^eV^ and 
I± = 0.86^^ = 3.28 X lO^^eV^, where the reduction of intensity for the _L beam through 
losses in the frequency doubling process has yet to be accounted for. 

((-L/ll)\ ^ 
Wq j TT is the area of the focal spot and the factor of 1/2 enters due 

to the splitting of POLARIS into two separate beams and uj±. It is clear, however, that 

these are upper theoretical estimates for the achievable intensities which will be certainly 

modified by the circumstances of the experimental setup. 

From these intensities, the peak electric field strength for the fundamental mode is -Ey = 

~ 1.43 X 10^ eV^. Further , it is appropriate to assume a relatively moderate conversion 
efficiency of 40% for the frequency doubled beam at these field strengths, yielding E± = 
VOATT w 3.62 X lO'^eV^. 

For high intensities, the external pulses as well as the probe beam must not only be 
spatially but also temporally well focused. Naturally, the pulses then have a spectral width 
Auj, which, for Gaussian pulses is related to the pulse length as Ao; ~ 0.44^. 

In order to detect the frequency shifted photons at tJout = '^in i with low noise, 
these photons should lie well outside the spectral widths Auj[^, Acjy and Alj_|_, centered 

extended to three spatial dimensions. Qualitatively, we expect that, for deviations from the orthogonal 
geometry, the resonance conditions Ak ~ and 5fc ~ receive an angular dependence based on the three- 
momentum transfer in all spatial directions. In consequence, the necessary condition uj± = 2tj|| is expected 
be modified as well. However, as this relation is experimentally easily accessible by frequency doubling, a 
beam geometry close to the orthogonal setup appears to remain viable. 
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around and and oj±, respectively. 

As the external pulses are comparatively long, we obtain a small spectral width of 
Awext — O.OleV. With the above pulse length, we see that Text ^ zf^^ as well as Text ^ Wq^^ 
are well obeyed if we assume similar focal properties for the fundamental beam mode and 
the frequency doubled beam. 

After the second upgrade, JETI++ is expected to provide an energy of around £ = 3J 
per shot at a central wavelength of around Am = 800nm (cjin = 1.55eV) with pulse lengths 
as small as Tin ~ 30fs. Assuming = 1, the smallest Rayleigh length for the probe beam 
is therefore ~ 12.7eV. From the pulse energy, the number of incoming photons per shot 
is given by Ni^ = £/(^m ~ 1-21 x 10^^. In addition, the requirements cJout ^ ^^\\ ' ^'^-L; ^'^H 
are well obeyed, since the spectral width is only Awjn ~ 0.06eV. 

As Tin ^ Text, all A'^in JETI photons are available for the conversion process as long as 
< zI,Wq, which is implemented above by the choice of the focusing geometry, see also 

Fig. g. 

A decisive experimental parameter is the pulse repetition rate which determines Aghot 
for a given measurement time. As the computation of the necessary statistics for the 
photon detection requires detailed knowledge of the laser specifications and the setup we 
perform our estimates for Aghot = Aout = 1- In particular, nonlinear processes within 
the experimental setup tend to modify the idealized Gaussian frequency spectra of the 
laser beams. If larger statistics for the outgoing photons are required, this can always 
be accommodated by a larger number of shots. In the present example, POLARIS, due 
to its higher energy, has the smaller repetition rate of both lasers, which is expected to 
approach /rep ^ O.lHz. In practice, O(IOO) shots per day can be achieved, being a huge 
accomplishment for a Petawatt-class laser. In the future, improved cooling schemes for the 
amplifying medium may even lead to a further enhancement of the repetition rate. 

For these parameters, the discovery potential follows from Eqs. (6.1-3.3). In Fig. |2|, 



we obtain two black wedge- like curves around the resonant masses "1-11= 2.73eV and m± = 
3.63eV, which are determined by the photon energies of POLARIS and JETI. The peaks of 
the wedges and thus the minimal accessible coupling strength lie at 5 ~ 7.7 x 10~^GeV^^ 
and g ^ 1.7 x 10~^GeV~^, respectively. 

The currently best laboratory limits on ALPs by the LSW setup of the ALPS collabo- 
ration [|lO|, are indicated as blue-shaded area in Fig. |2|, while the best limits on solar axions 
are provided by the CAST experiment [||, denoted by a green-dashed line. Although astro- 
physical considerations currently give the strongest constraints, they are somewhat model 
dependent due to a different momentum-transfer regime 

In dipole experiments, the external magnetic field B is essentially constant, and the 
axion-photon conversion and reconversion probability in vacuum are well approximated by 

/ gBL\^ sin^(x) 

for the ALP parameter space under consideration, cf. Fig. |2[ Above, we have abbreviated 
X = rin?L/{4:Uj^i^), where L denotes the length of the dipole magnet and oj^i^ is the axion 
and photon energy, respectively. From Eq. (|6.6|) it is obvious that the best exclusion bounds 
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are obtained for small arguments a; <C 1, i.e. for small axion masses at fixed energies Wt^y-y 
and dipole length L, cf. Fig. ||. If the axion masses become too large, the conversion 
probability suffers from a x'"^ suppression. 

As the CAST-experiment utilizes uj^ ~ keV solar axions, the best bounds are obtained 
for masses below m < leV. By contrast, in laboratory experiments generically uj-y ~ eV, 
and thus the drop-off sets in at even lower masses m ~ 10~^eV. In both situations, the 
accessible mass regions can be slightly extended to higher masses by the use of buffer gas. 
However, the generic form of the conversion probability for constant external fields, see Eq. 
(|6.6| ), disfavors dipole searches for the exploration of higher ALP mass ranges. 

For both dipole experiments, the exploration of the m ~ leV region is difficult and 
thus the purely laser-based search can complement the existing ALP searches in the large 
mass region as detailed above. An additional feature of the purely laser-based search is 
the strong sensitivity to the axion mass which originates from the condition of momentum 
conservation. On the one hand this can be advantageous, since it allows for a direct 
estimate of the axion mass if a signal is detected. On the other hand, it would of course be 
desirable to explore a larger range of the axion-mass-coupling plane within a single setup. 

For this purpose, it would be favorable if the involved laser frequencies were tunable 
within a certain frequency range. 

In fact this can be realized by the use of optical parametric amplifiers (OPAs), which 
are employed to tune the frequency of an optical pump laser over a wide frequency range 
while retaining the temporal structure of the pulse to a good approximation. However, for 
today's OPAs, the pump energy is limited to approximately ImJ. Since this requirement 
strongly limits the available intensity, we choose to tune the frequency of the probe beam 
rather than the external beam which enters twice in the setup in its fundamental and 
frequency doubled mode. For optical probe lasers, we have thus a limitation of the number 
of incoming photons to around A^in — 10^^. Yet, as A^in ~ 5^, cf. Eqs. (6.1-3.3), for the 
above considered setup the sensitivity to the coupling is reduced by just around one order 
of magnitude. 

A feasible tuning range for today's OPAs covers a large spectrum from infrared to the 
ultraviolet, see e.g. [^], which can in principle be further extended to larger frequencies by 
higher harmonic generation. Thus, asserting a tuning range of A < 2 x 10~^m {uj > 0.06eV), 
one could in principle explore axion mass ranges above m± > 2.46eV for cjout = Win-t-wy as 
indicated by the black solid vertical line in Fig. ^ for POLARIS intensities and frequency 
^. We find that, employing OPAs, the exclusion bounds on ALPs could be extended to 
g>1.8x 10-^GeV-^ 

However, it has to be taken care that the outgoing probe photons are shifted to fre- 
quencies outside the spectral widths of the employed lasers as discussed above. Also, the 
efficiency for the OPAs and the higher harmonic generation processes have in principle to 
be taken into account. In this order of magnitude estimate these effects are neglected. 

In order to estimate the discovery potential at future facilities, let us extend our con- 
siderations to the planned Exawatt facility ELI [pO|. At this facility, a potentially feasible 



®Here, we neglect the my solution, since for this mass the outgoing frequency is ujout = i^in — i^y and 
thus the frequency of the incoming probe photon cannot be reduced below uii . 
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intensity aim is at / = 10^^-^. According to the previous estimates for the available in- 
tensity, we thus obtain the ELI field strengths E\\ ~ 8.29 x 10^ eV^ and Ei_ ~ 5.25 x 10^ eV^. 
These field strength exceed those of the POLARIS/JETI setup by roughly two orders of 
magnitude. The achievable focal parameters can be expected to be of the same order of 
magnitude as in our previous estimate. 

In Fig. ^ the red dotted vertical line indicates the region which could be probed at ELI 
with an additional OPA with tuning range and A'^inA^shot — 10^^ as above. Here we find that 
the coupling region above (7 > 1.9 x 10~^GeV~^ could be explored. Intriguingly, already a 
moderate demand on the number of interacting laser photons can almost complement the 
bounds of ALPS in the higher ALP mass region. Ultimately, the red dot-dashed vertical 
line suggests the necessary requirements at ELI for tests of typical QCD axion models, 
which are plotted as a yellow band, see Q and references therein. In order to make contact 
with the range of QCD axion models, in the attainable mass region, A^inA^ghot ~ lO'^^ is 
needed at ELI to explore an ALP mass range up to g' > 3.4 x 10~^GeV~^. This is a rather 
strong requirement for current technology, but with the advance of OPA technology and 
in high-intensity laser technology this could be a worthwhile long-term aim for the future. 

7. Conclusion 

In conclusion, we have investigated the feasibility of a search for axion-like-particles in a 
purely laser-based setup. In particular, we have concentrated on modern high- intensity 
laser systems, since the available field strengths can serve as a lever arm for probing the 
weak coupling g of axions and axion-like-particles (ALPs) to electromagnetism. 

We have studied in a one-dimensional setting the conversion of a probe beam with 
frequency Win into an ALP beam in an external laser and its reconversion into a photon 
beam by a second external laser for Gaussian beam profiles in the formal limit of infinite 
pulse lengths. The important difference to conventional photon-axion mixing in homoge- 
neous fields arises from a split of the incoming photon into several partial waves whose 
frequencies are given by all non-negative sums and differences of the three laser frequency 
scales. 

We have also shown that the amplitude of these partial waves is proportional to two 
damping terms induced by the conservation of three-momentum at the conversion points. 
The constraints imposed by momentum conservation, together with the practical require- 
ment that the outgoing photon should be frequency shifted with respect to uj\a for reasons 
of detectability, demand for a carefully designed experimental setup. 

We have shown that these constraints can be satisfied in a specific set-up involving 
one external beam of frequency a;|| to counter-propagate with respect to the probe beam 
and another external beam of frequency ujx_ to propagate orthogonally to the probe beam. 
For the situation where uj± = 2oj\\, momentum is conserved at both conversion points at 
the same time while the frequency of the outgoing beam is different from the frequency 
of the incoming beam: ojout = + ^\\ or ojoxxt = '^in — ^\\ depending on the order of the 
interaction of the beams. 
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Figure 2: Axion-like-particle exclusion bounds for purely laser-based setups in comparison to 
searches with dipole magnets, see text for parameters and attainable coupling values. The blue- 
shaded area in the upper left corner gives the currently best laboratory bounds on axion-like- 
particles from the ALPS collaboration [|lO|, while the best limits on solar axions are provided by 
the CAST experiment Q , denoted by a green-dashed line. The black wedges denote the exclusion 
limits for a setup involving the JETI-|--|- and POLARIS laser systems at Jena for one shot at single 
photon detection. The black line indicates the principle exclusion bounds at this setup by frequency- 
tuning through optical parametric amplification. The red dotted line corresponds to an estimate 
for the best obtainable bounds with ELI with a present-day OPA system. The red dot-dashed line 
suggests the necessary requirements at ELI for testing the parameter regime of typical QCD axion 
models, which are plotted as a yellow band. 



This mechanism is reminiscent to sum-frequency generation and difference-frequency 
generation known from Nonlinear optics. 

The amplitudes of theses processes are peaked around certain resonant masses which 
are determined by momentum conservation: For the frequency up-converted probe photons 
with cjout = Win + the resonant mass is m_L = '^\J^'\\ + <^mW||, whereas my = 2^ui\ni^\\ 
for the down-conversion process with cjout = — ^\\ ■ 

We concluded that by frequency tuning by means of optical parametric amplification, 
purely laser-based experiments hold the prospect of providing the strongest laboratory 
bounds on axion-like-particles in the 0{eV) mass range. As laboratory searches with 
dipole magnets generically probe only lower mass ranges, purely laser-based experiments 
could complement them in an essential manner. 
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In summary, our work suggests that high-intensity lasers are about to evolve into a 
new tool for fundamental physics. 
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A. Appendix 



In the previous calculation, we have considered the interaction of two external lasers with 
frequencies uj± and with a probe beam Wjn, in which the first, orthogonally propagating 
beam mediates the photon-axion conversion and the second, counter-propagating beam 
mediates the back-conversion from axions into photons. 

However, if all beams are focused simultaneously onto the same spot as considered, 
it is of course not experimentally obvious, which of the beams causes the conversion and 
back-conversion process, respectively. For this reason, we want to consider the process at 
interchanged interaction order _L oH. 

To this end, we employ the beam parameterization of the counter-propagating external 
beam, cf. Eqs. (|5.2| ) and ( ^.3| ) for the photon- ALP-conversion process. Now, the electric as 
well as the magnetic component of the external field can interact yielding an extra factor 
of 2. Then, the ALP amplitude of Eq. ( [4.3| ) reads 



f'°° 1 1 

[z\t') = gE,^E\\ dz" = 

■'-^ Vi + (^74^)%/i+rz". 



dt".h {m^/{t' -t"Y -\z' - z"\'^'^ e {it' - t") - \z' - 
X sin ( Wint" — hmz" + arctan 



oo 



— oo 




We proceed analogously to the previous calculation and perform the temporal integration 



by the substitution t" — )• t' — T and by virtue of of Eq. (4.5). Specializing to |a;in±ti^||| > m 
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and using the identity in Eq. (4/7) one finds that the equivalent to Eq. ( [4.6D reads 



iz',t') = \gE;^E 



_5_g«('^in+<^||)<'g-isgn(2'-2;")fciz'g«V'|| 



°° 1 1 
X / dz" , „ I ■ ^ 7 



,i(-A:i„+fc||+sgn(z'-2")fei)z" 
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1 - q 



3i(-A:in-A:|i+sgn(z'-2")sgn(tJin-W||)fcax)2'' 



+ c.c. 



, (A.2) 



where here, according to the substitution _L -H-H, the axion wave vector becomes a function 
of the frequency of the counter-propagating external field: /c^ = {uj-m ± uj\\Y ~ "f^"^- 

Again, we find the characteristic structure of the ALP partial waves with frequencies 
■^ax = Win±W||, which have transmitted (sgn(z' — z") = +1) and refiected (sgn(z' — z") = — 1) 
contributions. The two remaining integrations over z" can be performed in the complex 
2;"-plane. The respective integrals read 



dz- 
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5k: 



6k: 



(A.4) 



where we have defined 



6k'^ = —km + + sgn(z' — z")k'^^ 

6k7 = —km — k\\ + sgn(z' — z")sgn(aJin — ^11)^3 



(A.5) 
(A.6) 



In Eqs. (A.3) and (A.4), we encounter again the resonant structure of the conversion 
amplitude which can be attributed to the requirement of momentum conservation. As 
before, the sharp momentum cutoffs as induced by the signum functions are in fact relaxed 
by an integration over a finite interaction region. In analogy to the previous considerations, 
we determine the resonant mass m. 
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In vacuum {uj = k), Eqs. ( A. 5 ) and ( [A.6D become 



-Win + Wll + sgn(z' - z")J (cJin + U}\\)'^ - ~ 



(A.7) 



(5fc|| = —Win — '^ii + sgn(z' — 2:")sgn(a;in — uj^^)y (win — W||)2 — ~ . (A. 8) 

The condition in Eq. ( [A.7D is solved by setting m = my = 2ya;inW|| in the case of trans- 
mission for Win > "^11 and reflection for Win < ^||- This can be understood intuitively: For 
momenta of the probe beam which are larger then the momenta of the counter-propagating 
beam, we find transmission, otherwise reflection of the ALP beam. 

As in the previously considered setup, for 6k'^ there exists no resonant mass. In the 
following, we thus again only keep the transmitted part of the sum- frequency solution w^. 



A justification for omitting the reflected part will be given below Eq. ( A. 14 ). 



We have as a pendant to Eq. (4.14): 



^m(z',t') « -^J-^5i?ini?||[ sin((win + W||)t' - ^^z' + 



1 - sgn(5fc+))e^''n' + (1 + sgn(,5A:+))e" 



(A.9) 



We now turn to the back-conversion of the ALPs into photons. However, as the external 
beam for the back-conversion propagates orthogonally to the z-axis, only the magnetic 



or electric field component can couple, cf. Eq. (5J.). However, due to the asymmetric 
coupling structure, Cout is not invariant under this choice, as discussed in Sec. |2[ 

As the ALP is massive, the contribution from the magnetic field component of the 
external beam will be larger, since it couples to the temporal derivative of (p. We thus 
choose = in the following. 

Following the steps below Eq. ( |5.1D , where is now given through Eq. ( |4.1| ), we find 
that after the integrations over t' and z' for the back-conversion, we end up with 



(1 - sgn{6k+))e^\\ + (1 + sgn(5A:+))e^ 



-5kt 



„i(win+W||+a;x)(i-sgn(z-z')^)p«(V'||+V'±)p-i(«'o"'5fc|)^ 



_ £gi(win+aj||-a;x)(t-sgn{2-2;')^)gi(V'||-'0±)g-i(«'o-'5fcx)^ + C.C. 



, (A.IO) 



with the abbreviations 



5k 
6k 



-k^^ + sgn(z - 2;")(win + Wy + W_l) 
-A;+ + sgn(z - z")(t^m + ^ii - w_l) . 



(A.ll) 
(A.12) 
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In Eq. ( [A. 10 ), we encounter the familiar behavior of the reconverted photons. The outgoing 



electromagnetic wave is composed of two partial waves whose amplitude is tied to the 
vanishing of the sum of the momenta in the process: 5k^. As before, we are interested in 
the situation of concurrent momentum conservation in both conversion processes, under 
the constraint Wout 7^ "^in- Thus, we rewrite Eqs. ( [A.ll ) and (A. 12) as 



—6k^ = (win + a;||)2 — — sgn(z — z'){LOin + ^\\ + ^ (A.13) 



I 



—5kj_ = (win + wy)^ — m? — sgn(z — z){ijj\n + 0J\\ — f^±) ^ , (A. 14) 

where we have multiplied the equations by —1, which is justified since appears as a 



square in Eq. ( A. 10 ) 



By comparing the above conditions for momentum conservation in the back-conversion 
process to the condition for the photon- ALP conversion in Eq. ( [A.7D , we again find, that for 
positive frequencies only JfcJ can be "matched" to the vanishing of 5k'^ for transmission and 
the relation uj_\_ = 2oj\\ . As the frequency of the outgoing beam evaluates to Wout = — ^\\ 
in this situation, we again find that only the transmitted parts of both waves can propagate 
undamped, since the solution requires ujin > (cf. discussion below Eq. ( |A.8| )). 

In summary, we find for the outgoing electromagnetic wave at uj± = 2a;|j the transmit- 
ted part 

(1 - sgn{6k))e^^'''~^ + (1 + sgn{6k))e~^'''^^] sin((a;in - a;||)(t - z) + V,, - V'x) , (A.15) 
where we have set Sk'^ = 6kJ^ = 5k. 
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